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ABSTRACT 

O 

^ . We present recent measurements of the orbital motion in the binaries, DF Tau 

p ^. and ZZ Tau, and in the triples, Elias 12, T Tau, and V853 Oph. We observed 

^ . these systems with the Fine Guidance Sensors on the Hubble Space Telescope 

and with adaptive optics imaging at the W. M. Keck and Gemini North Obser- 



> , 

2 ! T Tau Sa-Sb. Because the orbital coverage in most of these systems does not 

^ ! yet span a sufficient portion of the orbit, we are not able to find definitive orbit 

G\ 
O . 

\0 . parameter space allowed by the current set of data available for each binary. We 

^ ' constructed weighted distributions for the total mass of the binaries derived from 

a large sample of possible orbits that fit the data. These mass distributions show 
that the total mass is already well-defined. We compute total mass estimates of 



vatories. Based on our measurements and those presented in the literature, we 
perform preliminary orbital analyses for DF Tau, ZZ Tau, Elias 12 Na-Nb, and 



O 

H : O.78lg:?^M0, O.66l{]:JfM0, LUI^qUMq, and i.lStlfrMQ for DF Tau, ZZ Tau, 

j§ ! Elias 12 Na-Nb, and T Tau Sa-Sb respectively, using a distance of 140 pc. For 

Elias 12 Na-Nb, where the orbital coverage spans ~ 164°, we compute a prelim- 
^ I inary orbit solution with a period of ~ 9-12 years. By including an earlier lunar 

^ I occultation measurement, we also find a likely orbit solution for ZZ Tau, with a 

period of ~ 32 years. With additional measurements to continue mapping the 
orbits, the derived dynamical masses will be useful in constraining the theoretical 
tracks of pre-main sequence evolution. 

Subject headings: binaries: visual, stars: pre-main sequence, stars: fundamental 
parameters 
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1. Introduction 

Determining masses and ages of pre-main sequence (PMS) stars is important for un- 
derstanding the physical processes of star formation. For clusters of young stars, these 
parameters are commonly determined from a star's location relative to evolutionary tracks 
on a color-magnitude diagram. However, different sets of theoretical tracks yield discrepant 
results, particularly for low-mass (M < IMq) young stars, where the predicted masses and 
ages can differ by factors as large as 2-3 (Simon 2001, 2005). These discrepancies result 
mostly from differences in the choice of atmospheric opacities and the treatment of interior 
convection. 

Calibration of theoretical evolutionary tracks requires comparison with high precision 
dynamical mass measurements. Methods used to measure dynamical masses of pre-main 
sequence stars include: (1) observations of eclipsing double-lined spectroscopic binaries, 
where analysis of the radial velocity variations and the eclipse light curve yield the individual 
masses and radii of the component stars (Popper et al. 1987; Casey et al. 1998; Covino et 
al. 2000; Alencar et al. 2003; Stassun et al. 2004); (2) mapping the Keplerian rotation of 
circumstellar disks to determine the mass of the central star (Simon, Dutrey, & Guilloteau 
2000; Dutrey, Guilloteau, & Simon 2003); and (3) combining spectroscopic and spatially- 
resolved observations of binary star orbits to determine the component stellar masses and 
distance to the system (Steffen et al. 2001; Boden et al. 2005). 

Currently, there are only a small number of low-mass PMS stars accessible to each of 
these techniques. Observing eclipsing binaries, for instance, requires the orbit be inclined by 
~ 90° along our hne of sight, whereas single stars with rotating circumstellar disks must have 
disks large enough to detect and resolve with millimeter- wave interferometers. Moreover, 
visual and spectroscopic detections of young binaries tend to sample systems distributed 
in different regions of orbital parameter space. PMS spectroscopic binaries typically have 
periods measured on the order of a few years or less (Mathieu 1994), corresponding to 
radial velocities amplitudes on the order of tens of km s^^. Meanwhile, instrumental angular 
resolution limits the detection of visual binaries with similarly short periods. At the distance 
of the Taurus and Ophiuchus star forming regions (140 pc), the resolutions achieved through 
speckle interferometry and adaptive optics imaging at large 8-10 meter class telescopes or by 
the Hubble Space Telescope (HST) are capable of resolving visual binaries with periods on the 
order of decades. The current generation of long-baseline optical/IR interferometers brings 
visual binary measurements (Boden et al. 2005; Schaefer, Simon, & Prato 2005) further into 
the orbital parameter space of PMS spectroscopic binaries. 

Through a statistical comparison based on a sample of PMS dynamical mass measure- 
ments, only 9 of which are below 1 Mq, Hillenbrand & White (2004) showed that for masses 
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above ~ 1.2Mq, the evolutionary models reproduce the measured masses within la on av- 
erage. However, for masses below 1.2Mq, most of the evolutionary models systematically 
underpredict stellar masses by 10-30%, with the models of Baraffe et al. (1998) providing 
the most consistent results. The goal of the work presented in this paper is to increase the 
number of low-mass PMS stars with dynamical mass measurements, thereby strengthening 
the tests of the evolutionary tracks. 

In this paper, we present results from our ongoing program (Simon, Holfeltz, & Taff 
1996; Schacfcr ct al. 2003) to map visual orbits of low-mass PMS binaries using the Fine 
Guidance Sensors on HST and adaptive optics imaging at the W. M. Keck and the Gem- 
ini North Observatories. We present recent measurements of the separation of the young 
multiple systems DF Tau, ZZ Tau, Ehas 12, T Tau, and V853 Oph. For ZZ Tau and the 
close binary (Ehas 12 Na-Nb) in the Elias 12 triple, we present prehminary orbit solutions. 
For the remainder of the systems, it is not possible to determine a reliable set of orbital 
elements because our measurements do not yet span a significant fraction of the orbital 
period. However, even with the limited orbital coverage, useful mass estimates can be deter- 
mined based on the observed curvature in the measured orbit. In this paper, we refine the 
statistical technique developed in Schaefer et al. (2003) for determining the total mass of a 
visual binary when less than half of the orbit is mapped. We have improved the statistical 
analysis by incorporating a Monte Carlo search to sample the orbital parameter space and 
by weighting each orbital solution by its probability when analyzing the distributions 
of masses produced by the search. We present our analysis of the orbital solutions in § 4 
and discuss the dynamical masses computed from these results in § 5. We conclude with a 
comparison of our dynamical mass estimates with predictions from evolutionary tracks and 
a brief discussion on the variability observed in the multiple systems. 



2. Observations and Data Reduction 

The stellar properties of the young multiples observed are hsted in Table 1. The table 
specifies the target's HBC number, right ascension and declination, V-band magnitude, B-V 
color, K-band magnitude, spectral type, and luminosity. If a binary does not have a spectral 
type identification for each component, then the spectral type of the unresolved system and 
the total luminosity of the binary are given. The last column provides an ordered list of 
references for the values in the table. 
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2.1. Fine Guidance Sensor Observations 

The Fine Guidance Sensors (FGS; Nelan et al. 2004) are a set of three optical inter- 
ferometers on board the HST. Each FGS consists of a pair of Koesters prisms that detect 
the tih of the incoming wave front in two orthogonal (x and y) directions. For a given FGS, 
the normalized difference in counts emerging from each side of the prism, as a function of 
the tilt of the incoming wavefront, produces the interferometric fringe pattern (S-curve) of 
the target. For a binary system, the interferometric response is the linear superposition 
of two single-star S-curves, scaled by the relative brightness of the components and shifted 
according to their separation on the sky. 

We used FGSlr in its high angular resolution transfer mode to observe our targets in 
the F853W filter. Table 2 lists the dates of the FGS measurements presented in this paper. 
Each observation was obtained during one HST orbit and is composed of multiple (>10) 
interferometric scans. Scan lengths of 1.0" to 1.5" were used to measure the FGS baseline on 
either side of the binary. The raw x and y FGS scans were cross-correlated, co-added, and 
smoothed by piecewise polynomial fits to retain the interferometric signal but remove the 
effects of the statistical noise of counting photons (see Schaefer et al. 2003; Simon, Holfeltz, 
& Taff 1996). 

FGS observations of a binary can be deconvolved by modeling the system with the obser- 
vation of an unresolved single star. We selected the FGS single-star calibrators SAO 185689 
and HD 233877 because their colors (B-V = 1.50 and 1.10, respectively) closely approximate 
those of the targets. We used software distributed by the Space Telescope Science Insti- 
tute (STScI) to analyze the FGS observations of binary systems (Nelan & Makidon 2002). 
The routine follows a least-squares approach of searching through a grid of separations and 
brightness ratios for the two components (Lattanzi et al. 1992) and finding the minimum 
between the model and the FGS data. 

For triple systems (i.e. Ehas 12), we developed a program that follows a similar least- 
squares method. The routine searches over a grid of positions and brightness ratios for three 
stars. As with the STScI software for binaries, the best-fit solution is found from where 
the between the model and the FGS data is at a minimum. We searched down to step 
sizes of 0.1 mas in separation and 0.0001 in fractional flux. When applied to binary systems, 
the results derived from our program are consistent with those obtained with the software 
designed at the STScI. Figure 1 shows the model fit for the y-axis of the triple system, Elias 
12, on 2005 August 6. The left panel shows the residuals when Ehas 12 is fit as a binary 
(fitting the close pair as a single source), while the right panel shows the improvement when 
Elias 12 is fit as a triple. 
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The uncertainties in the positions and fractional fluxes are determined from the param- 
eter range that causes a variation of 1 from the minimum x^. The largest source of error is 
dominated by a mismatch between the S-curves of the target and calibrator. For the purpose 
of calculating the l-a uncertainty intervals, we scaled the so that the minimum is equal 
to the number of degrees of freedom, u (i.e. the reduced xl = X^/^^ is equal to 1). 

After the position of each component is measured from the FGS scans, the angular 
separation, p, of the stars on the sky is given by 

p = ^A9l + A91 (1) 

where AO^ and A9y are the projected separations on the x and y axes of the FGS. The 
position angle (PA) of the system, measured east of north, is determined from 

PA^2n-(f)+ arctan (-r^), (2) 

ABx 

where (p is the spacecraft roll angle (Holfeltz et al. 1995). 



2.2. Near-IR Adaptive Optics Observations 

We used the adaptive optics (AO) systems at the W. M. Keck and Gemini North 
Observatories to obtain high resolution images of the multiple systems. For the hierarchical 
triples (Elias 12, T Tau, V853 Oph), the wide single component is bright and lies within 
the isoplanatic patch of the close pair and therefore provides a simultaneous point-spread 
function (PSF) for measuring the separation and flux ratio of the close pair. Measuring 
the separation of the binaries requires using a nearby single star as a PSF reference. To 
minimize changes in the AO correction, we selected the single star DN Tau as a PSF because 
its brightness is similar to that of our targets. We observed the PSF star immediately before 
or after the science observations on the targets. 

At the W. M. Keck Observatory, we used the near- infrared camera, NIRC2 (Wizinowich 
et al. 2000), with the adaptive optics system on the 10-m Keck II Telescope. The NIRC2 
detector is a 1024x1024 Alladin-3 InSb array. The images were taken with the narrow-field 
camera which has a plate scale of 0.00994±0.0005" pixel"^ and an orientation of 0.7° ±0.2° in 
position angle (R. Campbell, private comm.). For each target, we obtained sets of dithered 
images, using a five-point dither pattern, with a dither offset ~ 2". 

We also obtained high resolution images with the near infrared imager, NIRI (Hodapp 
et al. 2003a,b), using the Altair adaptive optics system (Herriot et al. 2000; Veran et al. 
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2003) on the 8-meter Gemini North Telescope. NIRI contains a 1024x1024 AUadin InSb 
array. We used it in the f/32 configuration which has a plate scale of 0.0218" pixel"^ and 
an orientation of PA=0.0° ± 0.05° (Beck et al. 2004). For each target, we obtained sets of 
dithered images, using a five-point dither pattern, with a dither offset ~ 2". 

AO images of V853 Oph were also obtained with the sht-viewing camera, SCAM 
(McLean et al. 2000), during an observing run with the NIRSPEC spectrometer on the 
Keck II telescope. The SCAM detector is a 256 x 256 HgCdTe array. When used with the 
AO system, as we did, it has a plate scale of 0.017" pixel"^ (Wizinowich et al. 2003). 

Table 3 lists for each multiple system, the date of the AO observations, the instrument 
and filters used, AO correction rates, integration times, number of co-added exposures per 
image, and total number of images obtained. During the observations, the brightest com- 
ponent in each system was used as the AO guide star. The Strehl ratios, estimated by 
comparison with an Airy function, ranged from 0.1 to 0.3, depending on the seeing and the 
brightness of the target. The images were flatfielded using dark-subtracted, median filtered 
fiat fields obtained on the nights of observation. Pairs of dithered exposures of the objects 
were subtracted to remove the background. Figure 2 shows co-added AO images for T Tau, 
Elias 12, V853 Oph, DF Tau, and ZZ Tau. The close companions in Elias 12 and T Tau are 
easy to recognize. In the image of T Tau, T Tau Sa appears fainter than Sb (sec Beck et al. 
2004 and Duchene et al. 2005 for a detailed discussion of the photometric variability in the 
T Tau system). The close pair in V853 Oph A (the brighter component) does not appear to 
be resolved. 

For each of the multiple systems, we extracted subarrays from the image centered on 
the close binaries, with widths of 0.12"-0.30". For Elias 12, T Tau, and V853 Oph, we used 
the wide components (Elias 12S, T Tau N, and V853 Oph B, respectively) as the PSF. For 
DF Tau, and ZZ Tau we used separate observations of the single star, DN Tau, as a PSF 
reference. Since the NIRC2 observations of Elias 12 on 2004 December 24 were saturated on 
Elias 12N, we used an observation of DN Tau as the PSF for these images also. For all of the 
images, we used the PSFs to construct models of the close binaries by searching through a 
grid of separations and fiux ratios. We used the IDL INTERPOLATE procedure to shift the 
PSF by sub-pixel intervals using cubic convolution interpolation. To remove any excess light 
of the primary star from the extracted sub-array of the close pair, we followed a least-squares 
approach to fit the background as a slanting plane. We determined the position and flux 
ratio of the companions in the close binary from where the between the model and the 
observation reached a minimum. The uncertainties were determined by analyzing multiple 
images individually and computing the standard deviation. 
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3. Results 

Table 4 lists the results from the most recent FGS and AO measurements for DF Tau, 
ZZ Tau, Elias 12, T Tau, and V853 Oph. The separation and position angle, measured east of 
north, are given for each epoch of observation. The last column specifies the instrument used 
to obtain the measurement. For the triples, Elias 12 and T Tau, we report the separation 
measured for each of the wide pairings (Elias 12 S-Na, Ehas 12 S-Nb, T Tau N-Sa, and T 
Tau N-Sb), as well as the relative separation of the close pairs (Ehas 12 Na-Nb and T Tau 
Sa-Sb). Because of improvements in the method of data analysis and error estimation, we 
report updated results for the FGS observations of Elias 12 and the AO observations of Elias 
12 and T Tau that were presented previously in Schaefer et al. (2003) and Beck et al. (2004). 
We attempted to model the observations of V853 Oph as a triple system, however, it was 
not possible to obtain a reliable set of measurements for the close pair. As a result, we list 
in Table 4 only the relative position of the wide pair, V853 Oph A-B, fitting the close pair, 
V853 Oph A, as an unresolved source. 

The FGS records the number of counts detected by the photomultipliers (PMT) every 
25 msec. We used the flux ratios measured from the FGS scans to convert the average 
number of counts per 25 msec received from the target to the PMT counts contributed by 
each stellar component. These values are given in Table 5 for DF Tau, ZZ Tau, and Ehas 
12. Table 6 lists the flux ratios measured from the adaptive optics observations. For the 
triples, we list the flux ratio of each component relative to the widely separated member, as 
well as the flux ratio between the two components in the close pair. For the binaries, we list 
only the flux ratio between the two components. For each date of observation, we indicate 
the fllter used for the listed flux ratios. Technical speciflcations of the fllters are given in the 
instrument manuals at the observatory web sites. Since the flux of T Tau N has remained 
constant in recent years, with an average magnitude of 5.53 ± 0.03 at K (Beck et al. 2004), 
we convert the relative fluxes of T Tau Sa and Sb into K-band magnitudes and present these 
results in Table 7. 



4. Computation of Orbital Parameters 

Based on the cumulative orbital motion measured for each of the binaries, we calculated 
orbital flts by following the three-dimensional grid search procedure of Hartkopf, McAhster, 
& Franz (1989) and Mason, Douglass, & Hartkopf (1999). We explored the parameter 
space by stepping through ranges of the period, time of periastron passage, and eccentricity 
(P, T, e). The program locates the orbital solution for which the between the measured 
and calculated positions at the times of observation is at a minimum. 
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Input to the orbit fitting program includes a set of starting values {P,T,e), the corre- 
sponding step sizes and search ranges for each of the three parameters, and the measured 
separation and position angle at the times of observation. At each step through the three- 
dimensional grid, a least-squares fit is perfomed (Hartkopf, McAlister, & Franz 1989) to 
determine the Thiele-lnnes elements {A, B, F,G). The angular semi-major axis (a), incli- 
nation (i), position angle of the line of nodes and the angle between the node and 
periastron (lu) are computed from the Thiele-lnnes elements (i.e. Couteau 1981). The data 
are weighted by their respective measurement errors and the between the calculated 
{xc, Vc) and observed positions {xo, Vo) are computed from 



where the errors, and ay, are propagated from the measurement uncertainties listed in 
Table 4. The set of parameters that produces a minimum in the surface is retained as 
the best-fit orbital solution. 

Formal errors in the orbital elements can be derived using the general matrix method for 
multiple regression (Bevington & Robinson 1992). These error estimates involve linearizing 
the equations of orbital motion through a Taylor expansion terminated after the first order. 
For orbits that are poorly determined, higher order terms may become important. For such 
cases, it is necessary to examine the shape of the surface. In such cases, the precision 
of the computed orbital elements can be evaluated from the distribution of orbital elements 
that produces a variation of one from the minimum x^. 

For binaries where the orbital coverage spans less than 180°, it is not possible to deter- 
mine a reliable set of orbital parameters (Hartkopf, Mason, & Worley 2001). This is the case 
for most of the systems presented in this paper. In order to explore the x^ surface for each 
of the binary orbits, we performed a Monte Carlo search by modifying the three-dimensional 
grid search procedure so that the values of P, T, and e on each pass through the program are 
chosen at random. If a set of orbital parameters yields a value that is within Ax^ = 1 of 
the minimum, the solution is stored. After a specified number of possible solutions are found 
(N=10,000), the results are plotted and examined. The projection of the Ax^ = 1 surface 
onto each of the orbital parameters corresponds to the 1-a uncertainty interval, assuming 
that the parameters are independent (Press et al. 1992, Ch.l5). During this procedure, we 
scale the x^ values so that = 1 at the minimum. 

Because of the hmited orbital coverage for the binaries, the x^ surfaces turn out to be 
broad and flat in the region of the minimum. Table 8 lists the 1-a uncertainty intervals for 
the orbital parameters of DF Tau, ZZ Tau, Elias 12 Na-Nb, and T Tau Sa-Sb. Upper search 
ranges of P=500 years and e=0.99 were used when exploring the parameter space. If the 
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Ax^ = 1 interval extends beyond these search ranges, only the 1-a lower limit is quoted. 
The results for the individual systems are described in more detail below. In Section 5, we 
discuss how the broad range of allowable parameters influences estimates of the total system 
mass. 



4.1. DF Tau 

DF Tau is a classical T Tauri star (Herbig & Bell 1988) that was first resolved as a binary 
during a lunar occultation survey (Chen 1990). From angularly resolved HST STIS spectra, 
Hartigan & Kenyon (2003) identified the spectral types of the primary and secondary as 
M2.0 and M2.5, respectively. Since its discovery, the orbital motion of the DF Tau binary 
has been mapped through speckle interferometry (Ghez, Neugebauer, & Matthews 1993; 
Ghez et al. 1995; Thiebaut et al. 1995; Balega et al. 2002, 2004; Shakhovskoj et al. 2006), 
HST imaging (White & Ghez 2001), and FGS interferometery (Simon, Holfeltz, & Taff 1996; 
Schaefer et al. 2003). Preliminary orbit solutions were computed by Thiebaut et al. (1995), 
Tamazian et al. (2002), and Shakhovskoj et al. (2006), with periods ranging from 74 to 92 
years. However, the reliability of these results is limited by the sparse orbital coverage. As 
illustrated by Schaefer et al. (2003), a large range of orbital solutions can fit the available 
data. 

Figure 3 shows the orbital motion of DF Tau, using the cumulative astrometric data in 
the literature (referenced above) and our recent measurements presented in Table 4. The 
orbital coverage spans ~ 110°. By exploring the surface using the Monte Carlo grid search 
technique, we found a minimum located at P 95.8 years. However, since the 1-a interval 
extends from a period of 37 years to periods greater than 500 years, the precise orbital period 
is not definitive. Figure 3 shows three examples of possible orbital solutions for DF Tau at 
periods of 50, 75, and 100 years, which agree within Ax^ < 0.4. 

Figure 4 compares the current set of measurements for DF Tau to the orbital solutions 
computed by Thiebaut et al. (1995), Tamazian et al. (2002), and Shakhovskoj et al. (2006). 
We also plot the 95.8 year orbit that we found at the minimum of the surface through our 
Monte Carlo search. It is clear that the solutions of Thiebaut et al. (1995) and Tamazian et 
al. (2002) are no longer consistent with the current measurements. It is hkely that the orbit 
of Shakhovskoj et al. (2006) and our best-fitting solution will require substantial revision in 
the years to come. 

The orbital calculations for DF Tau indicate that the system has not been observed 
with sufficient coverage to determine a definitive set of orbital parameters. This is consis- 
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tent with its designation as a grade 4 orbit in The Sixth Catalog of Orbits of Visual Binary 
Stars (Hartkopf & Mason 2006). This grade assignment indicates that any orbital determi- 
nation should be considered preliminary, where substantial revision may be required for the 
individual elements. This agrees with the results presented here. 

4.2. ZZ Tau 

Observed as an unresolved system, ZZ Tau is classified as a classical T Tauri star with 
a spectral type of M3 (Herbig & Bell 1988). A lunar occultation observation (Simon et al. 
1995) revealed the binary nature of the system. The orbital motion of ZZ Tau has been 
monitored with the FGS since 1994 (Simon, Holfeltz, & Taff 1996; Schaefcr et al. 2003). A 
preliminary orbit analysis by Schaefer et al. (2003) indicated that a large range of orbits fit 
the available data. 

Figure 5 shows the orbital motion measured for ZZ Tau, based on the cumulative set 
of FGS and AO measurements (sec references above and Table 4). The orbital coverage 
spans ~ 120°. Three examples of possible orbital solutions, at periods of 40, 50, and 60 
years which agree within Ax^ < 0.3, are overplotted in the figure. The earliest two FGS 
separation measurements of ZZ Tau are the most discrepant from the orbital fits. In these 
observations, the FGS was oriented so that the binary remained unresolved along the x-axis 
and therefore the positions are derived from the combination of consecutive observations at 
different HST roll angles (see Schaefer et al. 2003). 

The Monte Carlo grid search technique reveals a nearly flat surface out to the search 
limit of P=500 yrs. The broad 1-a confidence interval reaches a shallow minimum at a 
period of ~ 180 years, however, a large concentration of solutions (30%) exist with periods 
between 30 to 100 years. The shallowness of the minimum is caused by limited coverage and 
sampling of the orbital measurements. It is therefore not possible to determine a definitive 
set of orbital elements for ZZ Tau based on the FGS and AO data alone. 

To improve the orbital coverage of ZZ Tau, we added the lunar occultation measurement 
observed by Simon et al. (1995) to our dataset. During the occultation of ZZ Tau in 1991.6, a 
projected separation of 29 mas was measured along a position angle of 244°. The orientation 
of the occultation is indicated in Figure 6. We used a standard Newton Raphson approach to 
compute a joint fit to the FGS, AO, and lunar occultation (LO) data. Table 9 lists the best 
fit orbital parameters, the of the fit, and the number of degrees of freedom, u. Figure 6 
shows a plot of the best fit orbit. Based on the 31.5 year period, the addition of the LO 
measurement increases the orbital coverage by ~ 100°. Continued FGS and AO observations 
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over the next two decades should provide complete coverage across the ZZ Tau orbit. 

4.3. Elias 12 

Ehas 12 (HBC 404, V807 Tau) is a hierarchical triple. The primary, Ehas 12S, is a 
classical T Tauri star (Herbig & Bell 1988) with a spectral type of K5 (Hartigan & Kenyon 
2003). The close pair, Elias 12N, is located ~ 300 mas to the north of Elias 12S and is 
currently separated by ~ 45 mas. Observed as an unresolved system, Ehas 12N, has a spectral 
type of M2.0 (Hartigan & Kenyon 2003). The orbital motion of the wide components, Elias 
12 S-N, with the close pair unresolved, was measured by Ghez, Neugebauer, & Matthews 
(1993), Leinert et al. (1993), and Simon, Holfeltz, & Taff (1996). Simon et al. (1995) resolved 
the projected separation of the close pair along the direction of a lunar occultation. Our 
FGS and adaptive optics observations to map the motion as a triple system began in 1999 
(Schaefer et al. 2003). 

Figure 7 shows the motion of Elias 12 Na and Nb as measured relative to the wide 
component, Elias 12 S. Figure 8 shows the motion observed for the close binary, where the 
position of Elias 12 Nb is plotted relative to Elias 12 Na. Over the 6 years of resolved FGS 
and AO observations, the close pair has moved ~ 164° relative to one another. Since these 
measurements span nearly half of the orbit, a preliminary orbit solution can be computed. 
Table 10 lists the best fit orbital parameters computed through a standard Newton Raphson 
method based on the FGS and AO data. The of the fit and number of degrees of freedom, 
V, are listed in the table. The orbit, with a period of 9.4 years, is overplotted in Figure 8. 
The orbital parameters and uncertainties from the formal fit are similar to those determined 
from the Monte Carlo search technique, which are presented in Table 8. 

To investigate the reliability of the orbit fit based on the FGS and AO data, we compared 
our results to the projected separation of 23 mas measured during the lunar occultation ob- 
servation along a position angle of 95° in 1992.2 (Simon et al. 1995). Adding this observation 
extends the orbital coverage to span more than one full cycle of the computed period. Be- 
cause the lunar occultation observation measured equal fluxes for Elias 12 Na and Nb, there 
is an ambiguity of 180° in the projected separation. The two possible orientations are shown 
as dashed hues in the panels of Figure 9. For each of these two possibhhties, we computed a 
joint fit to the FGS, AO, and lunar occultation (LO) data. The orbital parameters derived 
for each fit are listed in Tabic 10 and overplotted in Figure 9. While the unsealed for the 
P=8.65 year orbit is lower than that for the P=11.85 year by Ax^ = 0.65, the derived total 
mass of ~ 1.0 Mq for the 11.85 year orbit is more consistent with the M2.0 spectral type of 
Elias 12N than the ~ 2.2 Mq value for the shorter period. These values assume a distance 
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of 140 pc, the average distance to the Taurus star forming region (Kenyon, Dobrzycka, & 
Hartmann 1994). 

The three orbit solutions hsted in Table 10 show that there remains some uncertainty in 
the orbital parameters. However, with an orbital period of 9-12 years, a reliable orbit solution 
will be forthcoming over the next few years. We will continue to measure the orbital motion 
of Elias 12 in our ongoing observing programs. Additionally, since Elias 12 S provides a 
reference to measure the motion of the close pair about their center of mass, this could allow 
us to calculate the mass ratio of Elias 12 Na and Nb. Combined with the total mass from 
the visual orbit and an estimate of the distance, this would provide the individual masses of 
each component. 

4.4. T Tau 

T Tau, the prototypical young star (Joy 1945; Ambartsumian 1947, 1949), is a hier- 
archical triple located in the Taurus star forming region. The visible light primary, T Tau 
N, has a spectral type of KG (Basri & Batalha 1990). The infrared companion, T Tau S, 
located ~ 0.7" to the south of T Tau N, was discovered by Dyck, Simon, & Zuckerman 
(1982). Modeling of its spectral energy distribution suggests that T Tau S is the most lumi- 
nous component, having a mass of ~ 2.5Mq and a large visual extinction from circumstellar 
material (Koresko, Herbst, & Leinert 1997). T Tau S was found to have a close companion at 
a projected separation of ~ 0.05" through speckle interferometry (Koresko 2000). Spatially 
resolved spectra indicate an intermediate-to-early spectral type for the infrared companion, 
T Tau Sa, and a spectral type of Ml for the close companion, T Tau Sb (Duchene, Ghez, 
& McCabe 2002; Duchene et al. 2005). T Tau Sa and T Tau Sb remain optically invisible, 
down to a limiting magnitude of ~ 19.6 (Stapelfeldt et al. 1998). Silicate and water-ice 
features indicate high-levels of foreground material toward T Tau Sa and Sb (Ghez et al. 
1991; Beck et al. 2004). Additionally, CO absorption in the spatially resolved spectrum of 
T Tau Sa suggests the presence of a circumstellar disk viewed nearly edge-on and located 
within ~ 3 AU of the star (Duchene et al. 2005). Beck et al. (2004) discuss the cumulative 
measurements of the apparent motion of T Tau S relative to T Tau N. Since its discovery, 
the orbital motion of T Tau Sb relative to Sa has been monitored by Koresko (2000), Kohler, 
Kasper, & Herbst (2000), Duchene, Chez, & McCabe (2002), Furlan et al. (2003), Beck et 
al. (2004), Duchene et al. (2005), Mayama et al. (2006), and Duchene et al. (2006) ^ 



^We do not include the new astrometric measurements reported by Duchene et ah (2006) in our orbital 
analysis, but have added the reference for completeness. The new positions are consistent with the set of 
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Figure 10 shows the orbital motion of T Tau Sb relative to T Tau Sa, using the cumu- 
lative astrometric data in the literature (see references above) and our recent measurements 
listed in Table 4. We do not include the measurements reported by Mayama et al. (2006) 
because there appears to be a systematic offset of ~ 20 mas in the measured positions they 
present as compared to the rest of the astrometric data. Based on the set of infrared mea- 
surements, the orbital coverage spans ~ 80° in position angle. Overplotted in the figure 
are three examples of possible orbits, at periods of 20, 30, and 40 years, which agree within 
Ax^ < 0.1. The Monte Carlo grid search reveals that the 1-a surface is broad and flat 
and extends out to the search limit of P=500 years. The minimum in the surface occurs 
at a period of ~ 20.8 years. 

The northern and southern components of T Tau are both detected in radio emission. 
The southern radio source appears to be associated with the infrared position of T Tau Sb 
(Johnston et al. 2003, 2004; Loinard, Rodriguez, & Rodriguez 2003; Smith et al. 2003). By 
combining nearly 20 years of radio observations with the positions measured in the infrared, 
Johnston et al. (2004) fit an orbit, with a period of 25.5 ± 8.8 years, to the T Tau Sa-Sb 
system. The minimum that we find through the Monte Carlo search using only the 
infrared measurements is consistent with this value. Our results are also in agreement with 
the recent astrometric orbital analysis reported by Duchene et al. (2006). The high spatial 
resolution maps of Johnston et al. (2004) show that the southern radio source can actually 
be resolved into two components. However, the separation between the two radio sources is 
smaller than that measured for T Tau Sa and Sb in the infrared. The authors suggest that 
the radio emission could be due to magnetic reconnection processes in the environment of T 
Tau S, that is not coincident with either star. The position of the radio emission is useful 
for placing constraints on the orbital motion of the T Tau S binary, however, additional 
infrared and radio observations are necessary to further understand the connection between 
these measurements. 



4.5. V853 Oph 

As an unresolved system, V853 Oph is classified as a classical T Tauri star (Herbig & 
Bell 1988) with a spectral type of M2.5 (Bouvier & Appenzeller 1992). Chez, Neugebauer, & 
Matthews (1993) resolved the wide component of the triple system at a separation of ~ 400 
mas using speckle interferometry. In 1992, the primary, V853 Oph A, was discovered to 



measurements used in our analysis and because they span the same interval of time, would not significantly 
alter the results from our Monte Carlo search technique. 
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have a close companion at a separation of 13 mas along the projection of a lunar occultation 
(Simon et al. 1995), with a measured brightness ratio of 0.83 at K. Beginning in June 1999, 
we have monitored the system with the FGS and through AO imaging. 

The diffraction hmit of the Keck II Telescope at H is ~ 40 mas and that of HST at 
V is ~ 60 mas. This makes it difficult to measure rehably the small separation of the 
close binary in V853 Oph A. For the AO observations, we used the fainter tertiary, V853 
Oph B, as a simultaneous PSF reference. Since V853 Oph B is separated by only ~ 340 
mas from V853 Oph A and the system lies low in the southern sky at Mauna Kea, the 
observations require good natural seeing to obtain a well-resolved PSF, even with the use of 
AO. For example, during the N1RC2 observations on 2002 March 8 and 2004 June 15, it was 
not possible to extract the tertiary without including a significant amount of light from the 
primary. Resolving the close pair in V853 Oph with the FGS is also complicated by the small 
separation. For some FGS scans, there appears to be residual signal remaining after modeling 
the close pair as an unresolved source. However, measuring separations smaller than ~ 14 
mas on an individual axis of the FGS is not reliable, unless the magnitude difference can 
be constrained by a well-separated measurement on the other axis (Nelan et al. 2004). This 
was not the case for V853 Oph A. As a result of these complications, we were not able to 
derive a reliable set of measurements for the close binary in V853 Oph. 

In Figure 11, we show the motion of the primary, V853 Oph A, modeled as an unresolved 
source, relative to the wide tertiary, V853 Oph B. To help guide the eye, we plot a model orbit 
with a period of 600 years and a total mass of 1.0 Mq to approximate the motion of V853 Oph 
A relative to V853 Oph B over the course of the observations. The measurements are clearly 
scattered around the model orbit. This scatter could be caused by the orbital signature 
of the unresolved companion shifting the photocenter of V853 Oph A. We tried fitting a 
sinusoidal function to the residuals but were not able to derive a significant periodicity. With 
improvements in the sensitivity limits of ground-based infrared interferometers, resolving the 
close companion of V853 Oph A could prove to be an ideal target for interferometry. 

5. Total Mass Estimates 

The total mass of a binary is determined by a^/P^, with the semi-major axis on a 
physical scale. To examine the distribution of masses allowed by the current set of orbital 

measurements, we used our Monte Carlo grid search technique to explore the 3-0" confidence 
interval (Ax^ = 9) for DF Tau, ZZ Tau, Elias 12 Na-Nb, and T Tau Sa-Sb. As in the previous 
section, we searched for orbits out to periods of 500 years and eccentricities of 0.99. After 
finding 10,000 possible orbits that fit the data, we computed the total mass for each orbital 
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solution, Mfot — {a/lAUy / {lyr / Py Mq. To convert the angular separations to physical 
sizes, we used the average distance to the Taurus star forming region, d=140 pc (Kenyon, 
Dobrzycka, & Hartmann 1994). For T Tau, we used the VLBA parallax measurement of 
141.5 ± 2.8 pc (Loinard et al. 2005). The resulting mass distributions for DF Tau, ZZ Tau, 
Elias 12 Na-Nb, and T Tau Sa-Sb arc shown in Figure 12. In these histograms, each orbital 
solution is weighted by its probability. We truncated the histograms for DF Tau and T 
Tau at 4 and 10 Mq respectively; masses beyond these cut-off values are not representative 
of the measured luminosities and spectral types of the systems. This condition tends to 
eliminate highly eccentric, short period orbits from the sample and typically removes less 
than 5% of the 10,000 solutions in the 3-a interval. There is a clearly defined peak in each 
of the mass distributions. The last row of Tabic 8 lists the most probable value for the total 
masses of DF Tau, ZZ Tau, Elias 12 Na-Nb, and T Tau Sa-Sb. Since the mass distributions 
are asymmetric, the reported values represent the median of the distributions. The range of 
the quoted uncertanties include 34% of the values on each side of the median. For ZZ Tau, 
where we computed a preliminary orbit solution, the dynamical mass estimate based on the 
weighted Monte Carlo distribution agrees with the value derived from the formal fit within 
1(7. For Elias 12 Na-Nb, where we found two possible orbit solutions based on the combined 
AO, FGS, and LO data, the dynamical mass determined from the weighted Monte Carlo 
distribution favors the value derived from the 11.85 year orbit. 

To further examine the typical shape of the allowed parameter space. Figure 13 shows 
crosscuts through the surface for different sets of orbital parameters for DF Tau. The plot 
of angular semi-major axis versus the orbital period shows the strong correlation between a 
and P that produces the well-defined distributions for the total mass. The plots of the total 
mass versus e and P, as well as the plot of P versus e, show the long tail of short period, 
highly eccentric solutions that give rise to unphysically large values of the total mass. It 
is interesting to note that even if the chosen search ranges arc extended beyond P=500 
years, the calculated value for the total mass remains unchanged since the mass distribution 
fiattens out at large periods, converging to a value corresponding to the computed median 
mass. 

The result that poorly determined orbits can produce reliable values for the total mass 
was noted by Eggen (1967) and Hartkopf, Mason, & Worley (2001). The system mass is 
measured from the instantaneous radial acceleration and separation of the stars. Although 
the complete orbit may not be known, the orbital solutions must reproduce the projected 
radial acceleration of the system at the current separation, measured through the curvature 
of the apparent orbit. As a result, the value of (a^/P^) can be determined with much greater 
accuracy than the independent uncertainties in the orbital elements suggest. 
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6. Discussion 

6.1. Comments on Component Masses 

DF Tau is the only binary in our sample for which there is a spectral type and luminosity 
measured for both the primary and secondary components. Using the evolutionary tracks 
computed by D'Antona & MazziteUi (1997, DM97), Baraffe et al. (1998, BCAH98), Palla 
& Stabler (1999, PS99), and Siess, Dufour, & Forestini (2000, SDFOO) and the luminosities 
and spectral types listed in Table 1, we compute the predicted total mass of the DF Tau 
binary for each set of tracks. To convert the spectral types to effective temperatures, we use 
the temperature scale conversion of Kenyon & Hartmann (1995). The age of the DF Tau 
binary predicted from the different sets of evolutionary tracks ranges from 0.5 to 2 Myr. In 
Figure 14, we compare the masses predicted from the theoretical tracks with our dynamical 
mass estimate of 0.7S11q'^IMq measured for DF Tau. The error bars for the predicted 
masses are derived from an uncertainty of ±1 subclass in the spectral type determination. 
The predicted masses from each of the tracks agree with our dynamical estimate within 
~ 1.5(7. Based on the current uncertainties, it is premature to make a qualitative assessment 
of the evolutionary tracks based only on the total mass estimate of DF Tau. We are working 
on a program to determine the spectral types and luminosities of the close companions in 
ZZ Tau and Ehas 12 N. 

For an age of 1 Myr, the Ml spectral type of T Tau Sb (Duchenc, Ghez, & McCabe 
2002) corresponds to a predicted mass of 0.4-0.7 Mq, depending on the choice of theoretical 
tracks. Combined with the total dynamical mass measurement of A.IS'^qIjMq for the T Tau 
S binary, this range implies a predicted mass of 2.4-5.3 Mq for T Tau Sa. Thus, T Tau Sa 
may be a Herbig AeBe star caught at a very early phase of its evolution. With an estimated 
mass of 2.1 Mq for T Tau N (White & Ghez 2001), our dynamical results suggest that T Tau 
Sa is the most massive component in the system. This is consistent with the photometric 
and spectroscopic findings in the literature (Koresko, Herbst, & Leinert 1997; Beck et al. 
2004; Duchene et al. 2005) and the orbit fits based on radio and infrared data (Johnston et 
al. 2004; Duchene et al. 2006). 

6.2. Variability 

The two recent FGS flux measurements for DF Tau presented in Table 5 are consistent 
with the results in Schaefer et al. (2003), where we found that the flux of the DF Tau primary 
varies while the secondary remains constant. In Figures 15 and 16, we plot the changes in 
the FGS counts measured for the components in the Elias 12 triple and for V853 Oph A-B. 
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As with DF Tau, it appears that the primaries, Ehas 12 S and V853 Oph A, show the most 
variabihty in the system. It is interesting that in all three systems, the primaries display the 
activity associated with accretion from a disk. 

In Figure 17, we show the cumulative magnitude variations observed for the components 
of T Tau, based on our measurements presented in Table 7 and measurements reported in 

the literature (Koresko 2000; Kohler, Kaspcr, & Hcrbst 2000; Duchene, Ghez, & McCabe 
2002; Furlan et al. 2003; Beck et al. 2004; Duchene et al. 2005; Mayama et al. 2006). These 
results show that the large variability of T Tau Sa continues in the recent measurements, with 
variations of ~ 2 mag taking place over the course of a few months. Duchene et al. (2005) 
speculate that this variability may be caused by a clumpy, nearly edge-on circumstellar disk 
rotating around the central star. 



7. Summciry 

1. We report recent measurements of the orbital motion in DF Tau, ZZ Tau, Elias 12, T 
Tau, and V853 Oph using the Fine Guidance Sensors on HST and adaptive optics imaging 
at the Keck and Gemini Observatories. 

2. Based on the cumulative set of orbital measurements, we find a likely orbital solution 
with a period of ~ 11.9 years for the Elias 12 Na-Nb binary. Similarly, by combining our 
FGS and AO measurements of ZZ Tau with the published lunar occultation observation, we 
find a preferred orbital solution with a period of ~ 32 years. Our analyses for DF Tau and 
T Tau Sa-Sb show that we are not yet able to determine definitive orbital parameters for 
these systems. 

3. We developed a technique to determine dynamical mass estimates for binary orbits 
with limited orbital coverage. By using a Monte Carlo technique to search the orbital 
parameter space, we constructed weighted distributions for the total mass in DF Tau, ZZ 
Tau, Ehas 12 Na-Nb, and T Tau Sa-Sb. We find total masses of O.78l°:?^M0, O.66l°:J^M0, 
1.13lo;o9-^0) 4.13lo;97-^0, respectively, for these systems by adopting a distance of 
140 pc. 

4. The uncertainties in our dynamical mass estimates, combined with the limited knowl- 
edge of the luminosites and effective temperatures of the close companions preclude us from 
making a qualitative assessment of the theoretical tracks of PMS evolution based only on the 
masses presented herein. With continued orbit mapping and precise distance measurements 
obtained through astrometric missions such as SIM and GAIA, we expect to measure the 
total dynamical masses of the binaries to better than 5%. Determining the luminosity and 



-18- 



effective temperature of the the close companions will be the next step toward the goal of 
testing the evolutionary tracks. 
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Table 1. Summary of Target Properties 



Object 


HBC 


Coordinates (J2000) 


V 


B-V 


K 


SpT 


L(L0) 


Ref 


DF Tau A 


36 


04:27:02.80 +25:42:22.9 


12.43" 


0.75 


7.0 


M2.0 


0.56 


1,1,2,3,4 


DF Tau B 






13.10 


1.37 


7.9 


M2.5 


0.69 


1,1,2,3,4 


ZZ Taut 


46 


04:30:51.40 +24:42:22.2 


14.34 


1.52 


8.5 


M3 


0.7* 


5,5,6,5,7 


Elias 12 S 


404 


04:33:06.62 +24:09:55.2 


11.56 


1.18 


7.2 


K5 


2.3 


1,1,6,3,3 


Elias 12 Nt 






13.35 


1.24 


8.3 


M2.0 


1.04* 


1,1,6,3,3 


T Tau N 


35 


04:21:59.42 +19:32:06.4 


9.90 


1.19 


5.53 


KO 


12.5 


5,5,1,8,7 


T Tau Sa 










9.7" 






9 


T Tau Sb 










8.5" 


Ml 




9,10 


V853 Oph^ 


266 


16:28:45.29 -24:28:18.4 


13.3 


1.1 


8.0 


M2.5 




11,11,6,11 



^^Binary 

^Triple 

"Variable 

*Total luminosity of binary 

References. — (1) White & Ghez 2001; (2) Chen ct al. 1990; (3) Hartigan & Kcnyon 2003; (4) 
Hillenbrand & White 2004; (5) Herbig & Bell 1988; (6) Simon et al. 1995; (7) Strom et al. 1989; 
(8) Basri & Batalha 1990; (9) This work; (10) Duchene, Ghez, & McCabe 2002; (11) Bouvier & 
Appenzeller 1992 
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Table 2. Log of Fine Guidance Sensor Observations 



Object UT Date 

DF Tau 2003 Jan 23, 2003 Nov 19 
ZZ Tau 2003 Sep 21, 2005 Nov 06 
Elias 12 1999 Oct 03, 2000 Aug 27, 2000 Sep 04 

2000 Dec 31, 2001 Mar 04, 2002 Dec 22 
2002 Dec 26, 2003 Nov 08, 2003 Dec 23 
2005 Aug 06 

V853 Oph 1999 Jun 10, 2000 Apr 24, 2000 Jun 25 

2001 May 01, 2001 Jun 18, 2001 Jul 04 

2002 May 04, 2003 May 20, 2003 Jul 07 
2004 Jun 05 
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Table 3. Log of Adaptive Optics Observations 



Object 


Date 


Detector 


Filt 


AO Rate 


T 

mt 


No. 


No. 










(Hz) 


(s) 


Co-add 


Im. 


Elias 12 


2002 Mar 8 


NIRC2 


H 


272 


0.25 


10 


5 








K' 


272 


0.25 


10 


10 




2002 Oct 30 


NIRC2 


H 


140 


0.20 


20 


10 








LL 


1 90 


1 8 

U . -LO 


1 

-LU 


K 
o 








K' 


140 


0.50 


10 


10 








K' 


120 


0.18 


10 


5 




2004 Dec 24 


NIRC2 


H 


420 


0.18 


10 


10 








K' 


420 


0.18 


10 


10 




2005 Mar 1 


NIRI 


Hcont 


200 


1.00 


20 


25 








Kcont 


200 


1.00 


20 


25 




2005 Dec 6 


NIRC2 


H 


272 


0.18 


10 


5 








K' 


272 


0.18 


10 


5 


T Tau 


2002 Oct 30 


NIRC2 


K' 


154 


0.18 


10 


15 












n 30 


1 n 


1 5 




2003 Nov 18 


NIRI 


Kcont 


1000 


0.25 


10 


10 




2004 Dec 24 


NIRC2 


Kcont 


660 


0.18 


10 


20 




2005 Feb 9 


NIRI 


Kcont 


1000 


0.25 


10 


5 




2005 Mar 9 


NIRI 


Kcont 


1000 


0.20 


10 


5 




2005 Mar 24 


NIRI 


Kcont 


1000 


0.20 


90 


30 








H2(l-0) 


1000 


0.20 


90 


35 




2005 Oct 20 


NIRI 


Kcont 


1000 


0.20 


20 


10 




2005 Dec 6 


NIRC2 


K' 


90 


0.18 


10 


15 


V853 Oph 


2002 Mar 8 


NIRC2 


H 


79 


0.5 


10 


10 








K' 


79 


0.5 


10 


10 




2003 Apr 15 


NIRC2 


H 


74 


1.0 


10 


20 








K' 


74 


1.0 


10 


20 




2003 Apr 20 


SCAM 


K' 


91 


0.5 


5 


8 




2004 Apr 3 


NIRI 


Fell 


1000 


1.2 


10 


20 








Br7 


1000 


1.0 


10 


20 




2004 Apr 8 


NIRI 


Fell 


1000 


1.2 


10 


20 
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Table 3 — Continued 



Object 


Date 


Detector 


Filt 


AO Rate 




No. 


No. 










(Hz) 


(s) 


Co- add 


Im. 








Br7 


1000 


1.0 


10 


20 




2004 Jun 15 


NIRC2 


H 


55 


1.0 


10 


10 








K' 


55 


1.0 


10 


10 




2004 Jul 14 


NIRC2 


H 


55 


0.18 


10 


15 








K' 


55 


0.18 


10 


10 


DF Tau 


2004 Dec 24 


NIRC2 


Kcont 


235 


0.18 


10 


10 




2005 Dec 6 


NIRC2 


H 


161 


0.40 


10 


5 






NIRC2 


K' 


161 


0.40 


10 


15 






NIRC2 


K' 


84 


0.40 


10 


5 


ZZ Tau 


2004 Dec 24 


NIRC2 


K' 


73 


0.40 


10 


10 
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Table 4. Recent Measurements of PMS Stars in Multiple Systems 



Year 


p(mas) 


P.A.(°) 


Instrument 


DF Tau 


2003.061 


101.5±2.0 


258.6±1.2 


FGS 


2003.883 


105.1±2.0 


249.6±1.1 


FGS 


2004.981 


109.8±2.0 


247.0±1.1 


NIRC2 


2005.930 


110.9±1.6 


241.7±0.8 


NIRC2 


ZZ Tau 


2003.723 


62.8±1.5 


87.6±1.4 


FGS 


2004.980 


61.2±1.4 


74.9±1.3 


NIRC2 


2005.849 


61.7±1.5 


67.2±1.4 


FGS 


Ellas 12 S-Na 


1999.756 


282.9±2.1 


324.31±0.43 


FGS 


2000.655 


286.9±1.8 


323.56±0.36 


FGS 


2000.677 


284.9±2.7 


323.60±0.54 


FGS 


2000.999 


287.0±1.8 


322.83zb0.36 


FGS 


2001.170 


286.2±0.9 


322.12±0.18 


FGS 


2002.182 


285.4±1.2 


321.34±0.29 


NIRC2 


2002.829 


282.0±0.9 


319.86±0.25 


NIRC2 


2002.974 


281.4±1.5 


319.25±0.31 


FGS 


2002.985 


281.8±1.1 


319.60±0.22 


FGS 


2003.854 


271.6±1.0 


317.33±0.21 


FGS 


2003.977 


273.2±1.0 


316.96±0.21 


FGS 


2005.163 


268.9±1.0 


312.69±0.22 


NIRI 


2005.597 


260.5±0.7 


310.92±0.15 


FGS 


2005.930 


259.2±1.0 


310.31±0.27 


NIRC2 



Ellas 12 S-Nb 



1999.756 305.0±4.4 321.29±0.83 FGS 
2000.655 284.6±4.1 316.81±0.83 FGS 
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Table 4 — Continued 



Year p(mas) P.A.(°) Instrument 



2000. 


,677 


285. 


,2±4, 


.6 


317. 


,56±0. 


,92 


FGS 


2000. 


,999 


282. 


,9±4, 


.2 


315. 


,60±0. 


,85 


FGS 


2001. 


,170 


279. 


,4±3, 


.6 


314. 


,66±0, 


,74 


FGS 


2002. 


,182 


259. 


,5±1, 


,4 


314. 


,79±0, 


,34 


NIRC2 


2002. 


,829 


248. 


,5±0, 


,9 


313. 


,52±0, 


,26 


NIRC2 


2002. 


,974 


247. 


,0±2, 


.2 


312. 


,80±0, 


,51 


FGS 


2002. 


,985 


248. 


,9±2, 


.4 


313. 


,14±0, 


,55 


FGS 


2003. 


,854 


227. 


,5±2, 


.4 


312. 


,98±0, 


,60 


FGS 


2003. 


,977 


228. 


,8±2, 


.4 


313. 


,58±0, 


,60 


FGS 


2005. 


,163 


220. 


,8±1, 


,3 


311. 


,63±0, 


,33 


NIRI 


2005. 


,597 


214. 


,8±1, 


,3 


311. 


,32±0, 


,35 


FGS 


2005. 


,930 


214. 


,8±1, 


.0 


311. 


,56±0, 


,27 


NIRC2 



Elias 12 Na-Nb 



1999, 


,756 


27, 


,0±4, 


,9 


287, 


,8±10.4 


FGS 


2000. 


,655 


33, 


,7±4, 


,5 


226. 


,3±7.6 


FGS 


2000. 


,677 


30, 


,0±5, 


,3 


231. 


,2±10.2 


FGS 


2000. 


,999 


36, 


,2±4, 


,6 


222. 


,7±7.2 


FGS 


2001. 


,170 


37, 


,4±3, 


,7 


217. 


,8±5.7 


FGS 


2002. 


,182 


40, 


,4±1, 


,7 


188. 


,4±2.4 


NIRC2 


2002. 


,829 


44, 


,5±1, 


,1 


177. 


,9±1.3 


NIRC2 


2002. 


,974 


45, 


,4±2, 


,7 


177. 


,0±3.4 


FGS 


2002. 


,985 


44, 


,5±2, 


,6 


178. 


,6±3.4 


FGS 


2003. 


,854 


47, 


,9±2, 


,6 


158. 


,4±3.1 


FGS 


2003. 


,977 


46, 


,8±2, 


,6 


153, 


,8±3.2 


FGS 


2004. 


,979 


47, 


,7±1, 


,9 


139. 


,7±2.3 


NIRC2 


2005. 


,163 


48, 


,3±1, 


,6 


137. 


,6±1.9 


NIRI 


2005. 


,597 


45, 


,8±1, 


,5 


129. 


,0±1.9 


FGS 


2005. 


,930 


44, 


,7±1, 


,3 


124. 


,3±1.5 


NIRC2 
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Tablc 4 — Continued 



Year p(mas) P.A.(°) Instrument 



T Tau N-Sa 



2002.829 


697.6±2.1 


183.25±0.26 


NIRC2 


2003.881 


700.4±2.9 


183.15±0.55 


NIRI 


2004.979 


698.5±2.2 


184.74±0.27 


NIRC2 


2005.108 


704.0±1.6 


184.42±0.14 


NIRI 


2005.184 


702.8±2.4 


184.36±0.20 


NIRI 


2005.225 


703.7±0.9 


184.27±0.09 


NIRI 


2005.801 


703.3±1.2 


274.91±0.11 


NIRI 


2005.930 


703.0±5.7 


185.79±0.51 


NIRC2 


T Tau N-Sb 


2002.829 


684.8±1.6 


192.02zb0.24 


NIRC2 


2003.881 


673.87±0.89 


192.62±0.51 


NIRI 


2004.979 


662.13±1.9 


194.34±0.25 


NIRC2 


2005.108 


660.20±0.92 


193.89±0.09 


NIRI 


2005.184 


658.47±1.7 


193.79±0.16 


NIRI 


2005.225 


659.26±0.78 


193.80±0.09 


NIRI 


2005.801 


653.74±0.36 


284.50±0.06 


NIRI 


2005.930 


652.65±3.7 


195.19±0.38 


NIRC2 


T Tau Sa-Sb 


2002.829 


106.5±2.6 


284.5±1.4 


NIRC2 


2003.881 


116.5±3.1 


291.0±1.6 


NIRI 


2004.979 


119.6±2.8 


297.2±1.4 


NIRC2 


2005.108 


120.7±1.9 


300.36±0.88 


NIRI 


2005.184 


120.3±2.9 


300.63±1.40 


NIRI 


2005.225 


121.5±1.2 


300.42±0.58 


NIRI 


2005.801 


123.7ibl.3 


303.25±0.58 


NIRI 


2005.930 


121.9±6.8 


304.8±3.2 


NIRC2 



Table 4 — Continued 



Year p(mas) P.A.(°) Instrument 



V853 Oph A-B 



1999 


.440 


365. 


,98±1, 


.30 


95, 


.53±0. 


,20 


FGS 


2000, 


.313 


359. 


,12±0, 


.50 


95, 


.OlibO. 


,08 


FGS 


2000, 


.483 


358. 


.35±0, 


.50 


94, 


.80±0. 


.08 


FGS 


2001, 


.331 


353. 


,33±0, 


.63 


94, 


.71±0. 


,10 


FGS 


2001, 


.462 


351. 


.73±0, 


.63 


95, 


.22±0. 


.10 


FGS 


2001, 


.505 


351. 


,95±0, 


.50 


94, 


.90±0. 


,08 


FGS 


2002, 


.340 


347. 


,11±0, 


.67 


94, 


.69±0. 


,11 


FGS 


2003, 


.286 


344. 


,68±2, 


.00 


94, 


.79±0. 


,33 


NIRC2 


2003, 


.300 


351. 


,19±2, 


.57 


94, 


.04±0. 


,42 


SCAM 


2003, 


.382 


336. 


,97±0, 


.73 


93, 


.28±0. 


,12 


FGS 


2003, 


.514 


340. 


.06±0, 


.61 


93, 


.99±0. 


.10 


FGS 


2004, 


.256 


340. 


,15±0, 


.64 


92, 


.93±0. 


,11 


NIRI 


2004, 


.269 


340. 


,20±0, 


.50 


93, 


.OOibO. 


,08 


NIRI 


2004, 


.428 


333. 


,33±0, 


.51 


93, 


.81±0. 


,09 


FGS 


2004, 


.534 


337. 


,14±1, 


.40 


93, 


.72±0. 


,24 


NIRC2 
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Table 5. FGS PMT Counts 



Date 


PMT Counts (per 


25msec) ^ 


DF Tau 


A 


B 




2003/01/23 


1285.1 


668.2 




2003/11/19 


1018.7 


648.7 




// lau 


A 


J_) 




2003/09/21 


132.1 


43.8 




2005/11/06 


119.2 


43.7 




V853 Oph 


A 


B 




1 QQQ /fi /I n 


279.9 


74.7 




2000/4/24 


233.1 


60.7 




2000/6/25 


194.6 


79.1 




9001 /5/01 


185.6 


65.1 




2001/6/18 


229.7 


79.1 




2001/7/04 


193.9 


71.0 




2002/5/04 


260.3 


83.4 




2003/5/20 


287.0 


65.9 




2003/7/07 


291.8 


76.5 




2004/6/05 


207.0 


58.3 




Elias 12 


S 


Na 


Nb 


1 QQQ /I n 


1841.2 


261.3 


1 Q 


2000/08/27 


1848.8 


273.6 


144.2 


2000/09/04 


1939.1 


245.7 


117.7 


2000/12/31 


1899.1 


259.1 


87.8 


2001/03/04 


1937.0 


258.4 


81.6 


2002/12/22 


1877.1 


254.3 


166.7 


2002/12/26 


1585.7 


258.7 


153.7 


2003/11/08 


1970.0 


258.4 


120.7 


2003/12/23 


2043.0 


227.5 


127.1 


2005/08/06 


1904.8 


259.4 


135.6 
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^ Uncertainties are dominated by the 
square root of the number of counts derived 
from counting statistics. 



Table 6. Adaptive Optics Flux Ratios 



Date Filter Close Pairt Triple Flux Ratiot 



DF Tau 


"R / A 








2004/12/24 


Kcont 


0.475 








2005/12/06 


K' 


0.457 








2005/12/06 


H 


0.550 








ZZ Tau 




B/A 








2004/12/24 


K' 


672 








V853 Oph 


B/A 








2003/04/15 


H 


0.340 








2003/04/15 


K 


0.287 








2003/04/20 


K' 


0.270 








2004/04/03 


Br7 


0.227 








2004/04/03 


Fell 


0.275 








2004/04/08 


Br7 


0.197 








2004/04/08 


Fell 










2004/07/14 


H 


S54 








2004/07/14 


K' 


297 








Elias 12 


Nb/Na 


S 


: Na : 


Nb 


2002/03/08 


H 


0.720 


1.0 : 


0.388 


: 0.279 


2002/03/08 


K' 


0.692 


1.0 : 


0.399 


: 0.276 


2002/10/30 


H 


698 


1.0 : 


0.381 


• 266 


2002/10/30 


K' 


0.723 


1.0 : 


0.322 


: 0.233 


2004/12/24 


H 


784 


1.0 : 


0.710 


• 556 


2004/12/24 


K' 


0.748 


1.0 : 


0.606 


: 0.452 


2005/03/01 


Hcont 


0.689 


1.0 : 


0.393 


: 0.271 


2005/03/01 


Kcont 


0.731 


1.0 : 


0.308 


: 0.224 


2005/12/06 


H 


0.731 


1.0 : 


0.388 


: 0.284 


2005/12/06 


K' 


0.736 


1.0 : 


0.326 


: 0.240 


T Tau 




Nb/Na 


N : Sa : 


Sb 


2002/10/30 


Br7 


3.45 


1.0 : 


0.020 


: 0.069 


2002/10/30 


K' 


3.38 


1.0 : 


0.022 


: 0.074 


2003/11/18 


Kcont 


5.14 


1.0 : 


0.011 


: 0.054 
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Table 6 — Continued 



Date 


Filter 


Close Pair^ 


Triple Flux Ratio^ 


2004/12/24 


Kcont 


2.99 


1.0 : 


0.022 : 


0.065 


2005/02/09 


Kcont 


8.00 


1.0 : 


0.009 : 


0.071 


2005/03/09 


Kcont 


4.76 


1.0 : 


0.010 : 


0.049 


2005/03/24 


H2 


3.71 


1.0 : 


0.014 : 


0.052 


2005/03/24 


Kcont 


4.07 


1.0 : 


0.012 : 


0.050 


2005/10/20 


Kcont 


2.87 


1.0 : 


0.020 : 


0.058 


2005/12/06 


K' 


1.95 


1.0 : 


0.035 : 


0.067 



'^Typical uncertainties for the AO flux ratios are ±0.022 
for DF Tau B/A, ±0.047 for ZZ Tau B/A, ±0.006 for 
V853 Oph B/A, ±0.019 for Elias 12 Na/S, ±0.012 for 
Elias 12 Nb/S, ±0.001 for T Tau Sa/N, and ±0.002 for 
T Tau Sb/N. 
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Table 7. T Tau Component Magnitudes (assumes a magnitude of K=5.53 for T Tau N) 



Date 


Filter 


Sa 




Sb 




2002/10/30 


Br7 


9.78±0, 


.04 


8.44±0, 


.04 


2002/10/30 


K' 


9.67±0, 


.06 


8.35±0, 


.05 


2003/11/18 


Kcont 


10.47±0, 


.09 


8.69±0, 


.04 


2004/12/24 


Kcont 


9.69±0, 


.04 


8.50±0, 


.05 


2005/02/09 


Kcont 


10.66±0, 


.04 


8.40±0, 


.03 


2005/03/09 


Kcont 


10.50±0, 


.06 


8.81±0, 


.04 


2005/03/24 


H2 


10.16±0, 


.04 


8.73±0, 


.04 


2005/03/24 


Kcont 


10.30±0, 


.04 


8.78±0, 


.05 


2005/10/20 


Kcont 


9.76±0, 


.04 


8.62±0, 


.03 


2005/12/06 


K' 


9.18±0, 


.10 


8.47±0, 


.07 
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Table 8. Orbital parameters from the 1-a confidence intervals and dynamical masses from 

the weighted Monte Carlo distributions 



Parameter 


DF Tau 


ZZ Tau 


Ehas 12 Na-Nb 


T Tau Sa-Sb 


^(yr) 

T 
e 

a (mas) 
a;(°) 


>37 
1974-2022 
0.05-0.76 

>96 
128-176 


>31 

1994-1999 
0.35-0.88 

>60 
115-134 
111-129 
268-336 


q 4+1-2 

1991.16^°:^ 
47+0-10 

-0.11 

36.8+?-? 
147+?? 


>15 
1991-2001 
0.09-0.87 

>82 
11-60 


Mij^J' (Me) 






1 iq + 0.36 

-'-■-'-'J-0.09 


4 1Q+1-58 



Table 9. Orbital Parameters for ZZ Tau A-B 



Parameter 



AO-t-FGS+LO 



M( 



^(yr) 

T 

e 

a (mas) 



140pc ^ 



(Mo) 



31.5±5.3 

1994.26±0.85 
0.401±0.085 
61. ±10. 
133.9±7.6 
127.1±7.6 
270.±20. 
0.64±0.38 
20.62 
20 
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Table 10. Orbital Parameters for Elias 12 Na-Nb 



Parameter 


AO+FGS 


AO+FGS+LO 


AO+FGS+LO 


P (yr) 


9.4±1.0 


8.65±0.14 


11.85+0.12 


T 


1999.16±0.21 


1999.29±0.10 


1999.00+0.24 


e 


0.48±0.11 


0.570±0.031 


0.300+0.035 


a (mas) 


36.7±2.3 


39.1±1.0 


37.41+0.78 


in 


147.±14. 


137.2±2.8 


158.4+5.0 




97.±11. 


97.4±5.7 


-1.+27. 


o 


115. ±13. 


108.4±4.1 


39. +31. 


(Mo) 


1.53±0.45 


2.19±0.19 


1.026+0.068 




8.46 


8.78 


9.43 


V 


23 


24 


24 
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-400 -200 200 400 -400 -200 200 400 
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Fig. 1. — Left: Smoothed FGS scan of Elias 12, along the y-axis, on 2005 August 6. The 
dotted hne shows the best-fit binary model, using SAO 185689 as the PSF. The residual 
signature of the close companion can be seen at position of ~ —100 mas. Right: Elias 12, 
modeled as a triple system. The components in the close pair are separated by 30.4 mas 
along the y-axis of the FGS, with a flux ratio of 0.5. There is a significant improvement in 
the residuals of the triple star fit as compared to those of the binary model. 
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Toy N ^^^Bh'*' 

■ Jifk T Tou Sb 

T Tou So 

NIRI AO, K-cont 0.5 " 
24 Mor 2005 


Ellos 12 No 
* 

::,JiA Elios 12 Nb 

LliOS Iz •^ ^^^^H'; 

NIRC2 AO, H-band 0.2 " 
24 Dec 2004 


V853 Oph A 

V853 Oph B 

NIRC2 AO, K-bond 0.2 " 
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Fig. 2. — Co-added images of T Tau, Elias 12, V853 Oph, DF Tau, and ZZ Tau taken with 
adaptive optics imaging using NIRC2 on the Keck II Telescope and NIRI on the Gemini 
North Telescope. The camera, filter, and date of observation are listed in the panels for each 
image. 
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Fig. 3. — Orbital motion of DF Tau B (secondary) relative to DF Tau A (primary). The 
location of the primary star is marked by an asterisk. Examples of three possible orbital 
solutions for DF Tau at periods of 50, 75, and 100 years are overplotted. 
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Fig. 4. — Previously published orbits of DF Tau computed by Thiebaut et al. (1995, Th95 - 
red), Tamazian et al. (2002, Ta02 - blue), and Shakhovskoj et al. (2006, SH06 - green). We 
also include the orbital solution that we find at the minimum of the surface through our 
Monte Carlo search ("best fit" - black). The orbits are plotted at yearly intervals. 
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Fig. 5. — Orbital motion of ZZ Tau B (secondary) relative to ZZ Tau A (primary) as 
measured by our FGS and AO observations. The location of the primary star is marked by 
an asterisk. Examples of three possible orbital solutions for ZZ Tau, at periods of 40, 50, 
and 60 years, are overplotted. 



-44- 



100 



50 



-50 




CO 
D 












Q_ 
0) 
CO 

CO 

1 




2005.8 




2005.0 




"2003.7 ^ 







z 




- 2001.6 



100 



50 

E-W sep (mas) 



-50 



100 



Fig. 6. — Same as Figure 5, but includes the projected lunar occultation measurement of 
ZZ Tau from 1991.6 (Simon et al. 1995). The solid line shows the direction of the lunar 
occultation along a position angle of 244°. The dashed hne shows the projected separation 
of 29 mas measured along this direction. Overplottcd is the best fit orbit (P=31.5 yrs) based 
on the FGS, AO, and LO measurements. The position in the orbit at the time of the LO is 
marked by the shaded square. 
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Fig. 7. — Orbital motion in the Elias 12 hierarchical triple. The positions of Elias 12 Na 
(red symbols) and Nb (blue symbols) are plotted relative to Elias 12S (located at 0,0). 
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Fig. 8. — Orbital motion of Elias 12 Nb relative to Na (marked by an asterisk). Overplotted 
is the best fit orbit (P=9.4 yrs) based on the FGS and AO measurements alone. 
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Fig. 9. — Left: Same as Figure 8, but includes the projected lunar occultation measurement 
of Elias 12 from 1992.2 (Simon et al. 1995). The solid line shows the direction of the lunar 
occultation along a position angle of 95°. The dashed line shows the projected separation of 
23 mas measured along this direction. Overplotted is the best fit orbit (P=8.65 yrs) based 
on the FGS, AO, and LO measurements. The position in the orbit at the time of the LO 
is marked by the shaded square. Adding the LO observation extends the orbital coverage 
to span more than one cycle of the period. Right: Same as the left panel but showing the 
alternative possibility for the projected separation (rotated by 180°). Overplotted is the best 
fit orbit (P==11.85 yrs) based on the FGS, AO, and LO measurements. 
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Fig. 10. — Orbital motion of T Tau Sb relative to Sa (marked by an asterisk). Examples of 
three possible orbital solutions, at periods of 20, 30, and 40 years are overplotted. 
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Fig. 11. — Motion of V853 Oph A as measured relative to the tertiary, V853 Oph B (located 
at 0,0). The positions were determined by modeling the close binary, V853 Oph A, as an 
unresolved source. Overplotted is a model orbit with a period of 600 yrs and a total mass 
of 1.0 Mq. The large scatter of the measurements around the hypothetical orbit could be 
caused by the presence of the close companion of V853 Oph A. 
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Fig. 12. — Weighted distributions for the total mass of DF Tau (top left), ZZ Tau (top 
right), Elias 12 Na-Nb (bottom left), and T Tau Sa-Sb (bottom right). These results were 
obtained by searching for orbital solutions within the Sa confidence interval (A^^ = 9) by 
using our Monte Carlo grid search technique. The total mass (a'^/P^) is derived for each of 
the 10,000 orbits found within this interval. In the histograms, each solution is weighted by 
its probability. 
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Fig. 13. — Crosscuts through the surface for DF Tau. The color codes correspond to the 
Icr (red), 2a (blue), and 3cr (green) confidence intervals. 
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Fig. 14. — Total dynamical mass of the DF Tau binary compared with masses predicted 
from the evolutionary tracks of D'Antona & MazziteUi (1997, DM97), Baraffe et al. (1998, 
BCAH98), Palla & Stahler (1999, PS99), and Siess, Dufour, & Forestini (2000, SDFOO). The 
solid line shows where the predicted and dynamical mass estimates would be equal. 
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Fig. 15. — Variability of the FGS PMT counts measured for the three components in Ehas 
12. 
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Fig. 16. — Variability of the FGS PMT counts measured for V853 Oph A and B. 
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Fig. 17. — Cumulative K-band variability measured for T Tau Sa (closed circles) and T Tau 
Sb (open circles). We assume a constant magnitude of K=5.53 for T Tau N (dashed line; 
Beck et al. 2004). 



